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Abstract —The current formulation of the optimal scheduling 
of appliance energy consumption nses as optimization variables 
the vectors of appliances’ scheduled energy consnmption over 
eqnally-divided time slots of a day, which does not take into 
account the atomicity of appliances’ operations (i.e., the un- 
splittable natnre of appliances’ operations and resulting energy 
consnmption). In this paper, we provide a new formulation of 
atomic scheduling of energy consumption based on the optimal 
routing framework; the flow configurations of nsers over mnitiple 
paths between the common source and destination nodes of a 
ring network are used as optimization variables, which indicate 
the starting times of scheduled energy consumption, and optimal 
scheduling problems are now formulated in terms of the user flow 
configurations. Because the atomic optimal scheduling resnlts in 
a Boolean-convex problem for a convex objective function, we 
propose a successive convex relaxation technique for efficient 
calculation of an approximate solution, where we iteratively 
drop fractional-valued elements and apply convex relaxation 
to the resulting problem until we find a feasible suboptimal 
solution. Numerical results for the cost and peak-to-average ratio 
minimization problems demonstrate that the successive convex 
relaxation technique can provide solutions close to, often identical 
to, global optimal solutions. 

Index Terms —Atomic scheduling, convex relaxation, demand- 
side management, energy consumption schednling, optimal root¬ 
ing, smart grid. 


I. Introduction 

W E study the problem of scheduling electrical appliance 
energy consumption in residential smart grid. Our goal 
in revisiting this well-known problem of energy consumption 
scheduling (e.g., 0 -J9|) is to establish a new formulation 
of optimal scheduling problems where we take into account 
the atomicity of operations by household appliances, and 
to provide efficient solution techniques for the formulated 
optimal scheduling problems. By atomicity, we mean the 
unsplittable nature of appliances’ operations and resulting 
energy consumption. 

Note that the scheduling of electrical appliance energy 
consumption is a key to the autonomous demand-side manage¬ 
ment (DSM) for residential smart grid in optimizing energy 
production and consumption; the scheduling is based on smart 
meters installed at users’ premises and the two-way digital 
communications between a utility company and users through 
the smart meters, and the typical goals of DSM includes 
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consumption reducing and shifting, which lead into lower 
peak-to average ratio (PAR) and energy cost ID- 

Since the energy consumption scheduling was formulated as 
an optimization problem using energy consumption schedul¬ 
ing vectors as optimization variables representing appliances’ 
scheduled hourly energy consumption over a day Q, there 
have been published a number of papers on the subjects 
of appliance energy consumption scheduling and related 
billing/pricing mechanisms based on this formulation. For 
instance, the issue of optimality and fairness in autonomous 
DSM is studied in relation with billing mechanisms in Q, 
while the same issue is studied but in the context of user 
privacy in Q. The cost and PAR minimization problems, 
which are separately formulated in Q, are integrated into 
a PAR-constrained cost minimization problem in |j7). This 
integration is also further extended to take into account con¬ 
sumers’ preference on operation delay and power gap using 
multiple objective functions. In Q, instead of typical concave 
n-person games, Rubinstein-Stahl bargaining game model is 
used to capture the interaction between the supplier and the 
consumers through a retail price vector in lowering PAR to a 
certain desired value. 

In most existing works on the energy consumption schedul¬ 
ing for autonomous DSM in smart grid, however, no serious 
attention has been given to the microstructure of scheduled 
energy consumption over time slots; their major focus is on the 
optimal value of an objective function that is only based on the 
aggregate load from the scheduled energy consumption. Few 
exceptions in this regard include the works on the integration 
of consumers’ preference ||7| and the use of load consumption 
curves in the objective function 

Below are some scenarios illustrating the importance of the 
atomicity of appliance operations. 

• When a user washes clothes, the washing machine should 
be continuously on for a certain period depending on the 
amount of laundry, e.g., for two hours, during which the 
operation of the clothes washer cannot be interrupted and 
the supplied power cannot be reduced arbitrarily. As the 
washing task can be activated anytime within a specified 
period, e.g., from 9 AM to 3 PM when the user is out to 
work, the major goal of autonomous DSM is to determine 
the optimal two-hour time slot to complete the task when 
the energy price is lowest during the specified period. 

• For heavier-duty tasks like charging the battery of a plug¬ 
in hybrid electric vehicle (PHEV), again the continuity of 
the operation is important, e.g., four hours of uninterrupt¬ 
ible charging to maintain the lifetime of PHEV’s battery. 

• Appliances such as a rice cooker can significantly con- 
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tribute to the overall cost saving if handled properly. 
Some rice cookers’ function may take more than one 
hour for completion, e.g., the slow cooking function 
for delicate soup, whereby the cooking process can be 
scheduled within any time-shot during the day. 

Considering that most of the operations subject to DSM 
are either atomic as such (e.g., laundry cleaning by a clothes 
washer) or consist of atomic suboperations (e.g., house heating 
by a heater operating in the morning and in the evening) Q, 
therefore, we provide a new formulation of the optimization 
problem for atomic scheduling of appliance energy consump¬ 
tion and efficient solution techniques for resulting problems 
in this paper. Atomic scheduling has been mostly discussed 
in the context of concurrent task scheduling on a multi¬ 
processor/core system on a chip (SOC) (e.g., O)) or trans¬ 
action processing (e.g., ©)■ To the best of our knowledge, 
our work is the hrst attempt to formulate atomic scheduling 
of appliance energy consumption in the autonomous DSM for 
residential smart grid. 

The rest of the paper is organized as follows: In Section [H] 
we review and discuss the issues of the current formulation 
of appliance energy consumption scheduling based on energy 
consumption scheduling vectors dehned over time slots. In 
Section [III] we describe a new formulation of the appliance 
energy consumption scheduling based on the optimal rout¬ 
ing framework, which guarantees the atomicity of appliance 
operations, and a successive convex relaxation technique for 
the efficient solution of the Boolean-convex problem result¬ 
ing from the new formulation for a given convex objective 
function. In Section IV we demonstrate that the performance 
of successive convex relaxation technique through numerical 
results for the cost and PAR minimization problems. Section [V| 
concludes our work and discusses topics for further study. 


II. Review of Current Formulation of Appliance 
Energy Consumption Scheduling 
We hrst review the current formulation of appliance energy 
consumption scheduling problem by formally describing it. 
The formulation described here is largely based on |[T] but 
with some modihcations and extensions for clarity and better 
handling of scheduling intervals over a day boundary. Many 
notations and dehnitions in this section are applicable to the 
formulation of atomic scheduling problems in Section III as 
well. 

Let J\f= TV} denotes a set of users in a residential 

smart grid, whose appliances share a common energy source 
and subject to autonomous DSM. Without loss of generality 
and for ease of presentation, we assume that each user has 
only one appliance throughout the paperj^In this case a daily 
energy consumption scheduling vector of user n is dehned as 

( 1 ) 

where a scalar element denotes the energy consumption 
scheduled for a time slot h(^T-L= {0,..., iL—1}. A feasible 


^We use the term “user” and “appliance” interchangeably. 

^Time slot numbering in this paper starts from 0, which makes it easier 
handling energy consumption scheduling wrap-around the day boundary (i.e., 
from 11 PM to 6 AM) using modulo operation. See and for details. 


energy consumption scheduling set for user n is given by 



x^=E 


iT"<xl<iT\ xJj=o, yhenXHr. 


( 2 ) 


where 7 ™“ and 7 ™“^ are the minimum and the maximum 
energy levels for a time slot, En is the total daily energy 
consumption of user n’s appliance, and Tin is a scheduling 
interval dehned as follows 


Hn = {h\h = i mod H, ViG [a„, ^„]} (3) 

with a„G[0, H—1], 2H—2], and l</3„—a„<iL—1. 

With these dehnitions of scheduling vectors and feasible 
sets, the optimal scheduling is formulated as an optimization 
problem for a given objective function (e.g., total energy cost 
or PAR) as follows: 

minimize '0fx) (4) 

XnGVn, VtiGA/" 

where 

X = [xi,.. ., xn ] . (5) 

As shown in 0 . it is better to develop a DSM approach 
that optimizes the properties of the aggregate load of the users 
rather than individual user’s consumption. In fact, the optimal 
scheduling problems are formulated in terms of the total load 
across all users at each time slot h in ||T), i.e., 

Lh (x) ( 6 ) 

raGTV 

So Q can be expressed as 

minimize 6(L(x.)) (7) 

x„GA„. VnGV 

where 

i(x) = [Lo (x),... (x)]. ( 8 ) 


Note that the objective function becomes 

cl>{L{K)) = Y,Ch{LhiK)) (9) 

heV. 


for the energy cost minimization problem, where Ch{-) is a 
cost function indicating the cost of generating or distributing 
electricity energy by the energy source at a time slot h, and 


</>(i(x)) 


H maxL), (x) 
hen 

nGVV 


( 10 ) 


for the PAR minimization problem, respectively. 

A major issue with the formulation of optimal scheduling 
problems based on the energy consumption scheduling vectors 
in 0 is that it cannot guarantee the atomicity of appliance 
operations, which is illustrated in Fig. [T] (a). When a load 
from other appliances falls on the middle of the schedul¬ 
ing interval, especially that of non-shiftable appliances, the 
scheduled appliance energy consumption spreads over non¬ 
contiguous time slots, resulting in several gaps. Relatedly, the 
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Fig. 1. Examples of (a) non-atomic and (b) atomic scheduling of appliance 
energy consumption, where 7^“^, and 7 n^(-) are the minimum, the 

maximum, and the operating energy levels of an appliance for a time slot. 











6 i 2 3 4 5 ■■■ 


Fig. 2. An example of a predefined pattern of operating energy levels for 
Sn = 4. 


scheduled energy consumption may not provide enough power 
for appliances to carry out required operations because, during 
the scheduling, the energy levels over time slots are determined 
to achieve the optimal solution for the given objective function 
but not to meet the actual energy consumption requirements of 
the appliances for the operations. The atomic scheduling that 
will be described in Section on the other hand, assigns 
contiguous time slots with a predefined pattern of operating 
energy levels (i.e., 7^^(-)) as shown in Fig. [^(b), even at the 
expense of increased penalty in optimization. 

Ill. Atomic Scheduling of Appliance Energy 
Consumption 

We assume that all appliance operations are atomic (i.e., 
no gaps in energy consumption during the operations) and 
that the operation of the appliance of user n requires S„ 
contiguous time slots belonging to a scheduling interval of 
with a predefined pattern of operating energy level^ 7n^(h) 
such that 77(/i)>7;r*” Vhe[0, 5„-l] and 77 (/i)= 7 ™“ 
otherwise. Clearly, /3n>cin+<5n—1- Without loss of generality, 
we set 7 ™“ to 0; the resulting daily energy consumption En, 
therefore, is given by 

s„-i 

^—0 

Fig. shows an example of operating energy levels. 

Note that, if an appliance requires multiple, separate 
scheduling intervals (e.g., one for 6 AM-11 AM and the other 
for 1 PM-5 PM), it can be modeled as multiple (virtual) ap¬ 
pliances, each of them having only one contiguous scheduling 

^When the time slot duration is an hour (i.e., H=24), the energy level 
becomes a power level. 


interval (i.e., appliance 1 with 6 AM-11 AM and appliance 2 
with 1 PM-5 PM). 

A simple and straightforward formulation of atomic 

scheduling is using the starting times of operations as opti¬ 

mization variables, i.e., 

S= [si,...,SAr] (12) 

where a feasible set of starting times for user n is given by 

S^i ^ mod Ff, VzG[o;j.j, (13) 

The total load across all users at each time slot h, therefore, 
can be expressed in terms of starting times VriGAf, 

as follows 


L,(s) 4 ^77((/r-s„)modi/)/K„(,„)(/r) (14) 

n^J\f 

where 7?.„(s„) is a range of user n’s appliance operation for 
Sn defined by 

TZnisn) = {h\h = iuiod H, \/iG[sn, s„-|-(5„-l]} (15) 


and )(/i) is an indicator function; note that, for a set A, 

the indicator function is defined as 


lA{a) ^ 


1 if a G A, 
0 otherwise. 


(16) 


Because the feasible set is now discrete, we have to evaluate 
the objective function for all the elements in the feasible 
set, which makes the optimization problem impractical for 
large N and H. For instance, when A^=100 and H=2A with 
the worst case scenario of a„=0, I3n=‘23, and (5„=1 for 
all ugM, global optimization by direct enumeration would 
require evaluating the objective function 24^°° times, which 
is on the order of 10^^® times. To address the scalability issue, 
therefore, we provide an alternative formulation of atomic 
scheduling based on the optimal routing framework, which is 
amenable to convex relaxation and enables systematic analysis 
of suboptimal solutions with upper and lower bounds. 


A. Optimal Routing-Based Formulation 

We provide a new formulation of the optimal scheduling 
of appliance energy consumption based on the framework of 


optimal routing in networking 113 Ch. 5]; this new formula¬ 
tion takes into account the atomicity of appliance operations in 
scheduling without taking any additional objective functions. 

As for optimization variables, instead of the energy con¬ 
sumption scheduling vector x„ in ([T]), we use flow conflgura- 
tions of users over multiple paths between the common source 
and destination nodes of the network shown in Fig. 
which are defined as follows: 

The users belonging to the set J\f share a set V of paths 
connecting the source node S and the destination node D, 
where the path set V is defined as 


V 4 {p*-: 




(17) 


where i and j are the starting node number on the ring of 
a path and the number of hops, respectively; for instance, as 


^This is for the case of hourly time slots, i.e., H=24. 
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with starting time slot s and duration of slots. Fig. (b) 
shows mapping of all possible atomic appliance operations 
of two users (i.e., q;i=0, /3i=5, Si=2 for user 1 and a2=9, 
/32=14, (52=3 for user 2) into corresponding groups of flows 
(i-e., /i, • ■ ft for user 1 and /|,..., for user 2). 

We use the flow configurations of all users as the optimiza¬ 
tion variables, which are defined as follows; 

f 4 (18) 


where 

(19) 


A feasible atomic energy consumption scheduling set for user 
n is given by 


T =J { 


H /;=i. 


S^Sn 


( 20 ) 


/^e {0, 1}, Vse5„, /*=0, VseH\5„ 


where is the feasible set of starting times for user n that 
is already defined in ( [T3] l but in terms of s instead of Sn- 
Note that the constraint of ^ /*=! in @ ensures that 

an optimal user flow configuration vector /„ has only one 
non-zero element representing the starting time of a scheduled 
atomic operation. 

With the user flow configurations, the total load across all 
users at each time slot hGH can be calculated as the sum of 
all flows passing over the link mod h)^ j 

U{{) ^ Y. iih-s) modi?) ( 5] rjR^(s){h) 

\s£Sn 

( 21 ) 

Now that we have an expression of the total load at each time 
slot in terms of flows representing atomic operations of appli¬ 
ances, we can formulate the problem of optimal scheduling of 
appliance energy consumption for a given objective function 
of (/) {L (f)) with L (f) = [Lq (f)) ■ • ■) Lh -1 (f)] as follows: 

minimize <h(L{i)). (22) 

f„GJ=-„.VnGA^ rv w; 



Fig. 3. Atomic energy consumption scheduling based on optimal routing: 
(a) A network connecting the source (S) and the destination (D) through 
24 intermediate nodes with a sample path and its constituent links 

and (b) mapping of all possible atomic operations of 

two appliances into two groups of flows (/^,..., ff and /I,. •., over 
multiple paths on the network. 


illustrated in Fig. ^ (a), is a three-hop path consisting 
of links ^^d (^^^2 qjj jj^g j.jjjg ]s,fQjg ^g 

not take into account the radial links connecting the source 
and the destination nodes to and from the intermediate nodes 
on the ring in calculating hop counts and path/link costs 
We now define /® as a flow that user n sends on path 
which represents an atomic operation of user n’s appliance 


For a convex objective function, this problem becomes a 
Boolean-convex problem, since the sum constraint in ( |20l i is 
linear and the optimization variable /® is restricted to take 
only 0 or 1 as in IB- Note that this problem still cannot be 
solved efficiently but, unlike the problem formulation directly 
based on starting times in ([T^, is amenable to the successive 


convex relaxation technique that we describe in Section 111-B 


B. Successive Convex Relaxation 

We can obtain the convex relaxation of the atomic opti¬ 
mal scheduling problem in ( |2^ by replacing the nonconvex 
constraints /®e{0, 1} in ( |20| ) with the convex constraints 
as follows: 


^This means that their link costs are zero. 


minimize (j){L{i)) 

t^e:F„yneAf 


( 23 ) 
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E 


SGSn 


(24) 


0 < /^ < 1, Vse5„, /^=o, VseH\5„ L 


For a convex objective function, this problem becomes 
convex because the new feasibility set is now convex 
where all the equality and inequality constraints on f are linear. 
This convex-relaxed problem in ( |2^ , therefore, can be solved 
efficiently, for instance, using the well-known interior-point 
method 05). 

Let f denote a solution to the relaxed problem in ( |2^ . The 
relaxed atomic scheduling problem in ( [2^ is not equivalent 
to the original problem in because the elements of the 
optimal solution f can take fractional values (e.g., 0.75). 
The optimal objective value of the relaxed atomic scheduling 
problem in ( |2^ , however, provides a lower bound on the 
optimal objective value of the original problem in ( |22l l; the 
optimal value of the relaxed problem is less than or equal 
to that of the original problem because the feasible set for 
the relaxed problem contains the feasible set for the original 
problem. 

We can also use the solution of the relaxed problem in 
( |2^ to generate a suboptimal solution. Note that we cannot 
simply choose N largest elements of the optimization vector 
as is done for the optimal selection of sensor measurements in 
0; because the optimization variable f is a vector of vectors 
(i.e., /„), the selection of largest elements should be done per 
component vector. Also, we found that the convex relaxation 
spreads the starting times of appliance energy consumption 
over the same scheduling interval when there are several 
appliances with identical parameter values (i.e., a„, /3„, 
and 5n), which often results in many elements with the same 
values in the component vectors. If we choose randomly one 
element among the same-valued ones to break ties in this case, 
the performance of resulting suboptimal solution is not close 
to the optimal one. 

Therefore, here we present a new technique, called succes¬ 
sive convex relaxation, where we iteratively drop fractional¬ 
valued elements and apply convex relaxation to the resulting 
problem until we find a feasible solution. This technique 
is similar to the cutting-plane algorithm in mixed integer 


linear programming 116 Ch. 9] in that, at each iterative step, 
new constraints are added to refine the feasible region. The 
proposed technique, however, is much simpler in forming 
new constraints where it does not introduce any new slack 
variables and takes into account the structure of the opti¬ 
mization variables. The detailed procedure is described in 
Algorithmic where we introduce two variables, i.e., 6 *d, the 
threshold value for dropping, and Njj, the maximum number 
of fractional-valued elements that can be dropped per iteration. 
Note that the overall complexity of the interior-point method to 
solve the relaxed convex optimization problem is O 
operations 1151, which is needed per iteration and a dominating 
factor of the complexity of the proposed successive convex 
relaxation technique. With these variables, therefore, we can 


1 D 0 /* initialize the set of elements to drop */ 


2 

while true 


3 

f Solution of l|23^ with the equality constraints replaced 


by “fn = 0, Vs G 

C 

4 

Exclude elements that are 

the maximum of Vn G N, 


arrange in ascending order the remaining elements of f 
whose values are less than one, and let A, 12 ,... denote 


their indexes 


5 

fii -<—0 /* always drop the minimum-valued element */ 

6 

V 


7 

for j = 2,..., Nd 


8 

if hj < Od 


9 

fi, ^ 0 

/* drop it */ 

10 

V 


11 

else 


12 

Break 

/* return to the while loop */ 

13 

if there is only one nonzero element in / , VnGA/” 

14 

Break 

/* stop here; solution found V 


Algorithm 1: Successive convex relaxation. 


do a fine control of the number of fractional-valued elements 
dropped per iteration and the number of iterations to finish the 
procedure: A reasonable value of Qy) (e.g., 0.1) can prevent 
many high fractional-valued elements (e.g., >0.5) from being 
dropped unnecessarily per iteration when Nd is set to a rather 
large number. On the other hand, when a relaxed convex 
optimization problem gives a solution with a large number of 
fractional-valued elements smaller than Oo, we can drop them 
up to iVi). In this way we can speed up the whole procedure 
while not sacrificing the quality of the resulting suboptimal 
solution. 

C. Examples 

Below we provide specific examples of the optimal atomic 
scheduling for frequently used objective functions for DSM in 
residential smart grid. 

1) Energy Cost Minimization: For energy cost minimiza¬ 
tion, we can formulate it as the following optimization problem 
similar to 0: 

minimize Ch (Lhii)) (25) 

L6>'„,VnGA^ ^ 

hen 

where Ch{-) is a cost function indicating the cost of generating 
or distributing electricity energy by the energy source at a time 
slot h. 

Replacing the feasible set with in ( [25l ), we obtain 
the relaxed energy cost minimization problem as follows: 

minimize Ch (Lhii)) ■ (26) 

f„G.F„,VnGA^ ^ 

2) Peak-to-Average Ratio Minimization: Considering the 
total energy consumption J^neJV fixed, we can formulate 
the PAR minimization problem as follows: 

minimize max(L;,(f)). (27) 

Note that 0 is difficult to directly solve due to the max(-) 
term in the objective function. As mentioned in 0, however, 
this optimization problem can be turned into a Boolean-linear 
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TABLE I 

Appliance energy consumption requirements (adapted from 0) 


Appliance 

Parameters 

a [h] 

d [h] 

7°P [kWh] 

i [h] 

Dish Washer 

0 

23 

0.7200 

2 

Washing Machine (Energy-Star) 

0 

23 

0.4967 

3 

Washing Machine (Regular) 

0 

23 

0.6467 

3 

Clothes Dryer 

0 

23 

0.6250 

4 

Plug-in Hybrid Electric Vehicle 
(PHEV) 

22* 

29* 

3.3000 

3 


* Scheduling interval of 10 PM-5 AM. 


program, a special case of Boolean-convex optimization, by 
introducing a new auxiliary variable F as follows: 


minimize 


r 


subject to r > Lh{i), V/i S 'H. 


(28) 


Again, replacing the feasible set with in ( |28l ), we 
obtain the relaxed PAR minimization problem, i.e., 

minimize F 

r,f„G/'„,VnGAf (29) 

subject to F > Wh € %. 

Like the convex optimization, the linear program can be 
efficiently solved by either the simplex method or the interior- 
point method lig. Note that, when applying the successive 
convex relaxation technique described in Algorithm [T] the 
auxiliary variable F is not subject to dropping unlike f. 


IV. Numerical Results 


We demonstrate the application of the successive convex 
relaxation technique to the atomic scheduling of appliance 
energy consumption through numerical examples for the min¬ 
imization of total energy cost and PAR in the aggregated 
load. In all the examples, a day is divided into 24 time slots 
(i.e., H=24), and each user is to have an appliance ran¬ 
domly selected from the appliances whose energy consumption 
requirements are summarized in Table |Ij for simplicity, we 
assume constant operating energy levels for all appliances. 

For the energy cost minimization problem described in 
Section III-Cl we assume a simple quadratic hourly cost 
function as in Q, which is given by 

Ch {Lh) = anLl [cent] (30) 


where 


r 0.2 

I 0.3 


The objective function in ( |26] l 
be expressed as 


if/iG[0, 7], 

if/iG[8, 23]. ^ 

for this hourly cost function can 


^{i) = cj,{L{i))=Y,au{Luii)f. (32) 

h&'H 


Then its first and second derivatives, which are needed for 
convex optimization, are given as follows: For n€A7 and sSTf, 
the first derivative (i.e., the gradient of i/'(f)) is given by 


dm ^ 

dfn 

2 E {{h-s) modiF) (f) /k„(,)(L). 
h^n 


(33) 




(b) 


Fig. 4. Comparison of upper (UB) and lower bounds (LB) with true values 
from global optimization (GO) in atomic energy consumption scheduling for 
(a) cost minimization and (b) PAR minimization. 


For ni,n 2 (zAf and si,S 2 &'H, the second derivative (i.e., the 
Hessian of which is positive due to the convexity of 

the objective function in ( |32| ), is given by 

9^i/>(f) 

2 modiF) 7 °^ ((/i-Sz) modi?) (34) 

/tG« 

To evaluate the performance of the successive convex re¬ 
laxation technique, we first obtain the lower bound {LB) 
using ( [26l ) and ( |29l l for the minimization of energy cost 
and PAR, respectively. Then we obtain suboptimal solutions, 
which are also upper bounds, using Algorithm [T] with the 
dropping threshold 6d fixed to 0.1 and different values of No 
( UB{No))- We compare the lower and the upper bounds with 
the optimal objective value from global optimization {GO) 
using direct enumeration. The results are shown in Fig. 
where, due to the huge size of the feasible set for direct 
enumeration given in ( |20| l, we limit the maximum value of N 
to 10; for 7V=10, the size of feasible set is about 2.3x10^^, 
and in case of the energy cost minimization, it took 105.2 
hours (i.e., more than 4 days) to obtain the optimal solution 
from direct enumeration using an OpenMP-based parallelized 











































































































7 


version of C++ program on a workstation with two Intel® 
Xeon® processors running at 2.3 GHz providing 20 cores 
and 40 threads in total, while it took 23.3 seconds to obtain 
the suboptimal solution from the successive convex relaxation 
with 7V£)=1 (requiring 207 iterative steps of convex relaxation) 
using a MATLAB® script with OPTI toolbox on the same 
machine. Note that due to the random selection of appliances 
and their different requirements for energy consumption, the 
resulting energy cost and PAR are not proportional to N. 

The results in Fig. show that both lower and upper bounds 
are very close to the true optimal values. In case of the upper 
bounds, they are even identical to the true optimal values. 
For instance, the upper bounds on the energy cost for all 
values of Njy when N=2 and for N£)=l,2,5 when N=5 
are identical to true optimal values; in case of PAR, the 
upper bounds for all values of No when N<6 are identical 
to true optimal values. Considering the huge difference in 
the computational complexity between the two approaches, 
these results are remarkable. Of the results for energy cost in 
Fig. 1^ (a) and PAR in Fig. ffl (b), we found that the integrality 
gap in convex relaxation |1^ is more visible for the lower 
bounds on the PAR: Note that, unlike the total cost whose 
calculation involves all the appliances, the PAR depends on 
less number of appliances in its calculation, i.e., only those that 
contribute to a specific time slot whose load is maximum. The 
relaxation process, however, makes more appliances contribute 
to the PAR calculation by spreading fractional-valued flows 
over all possible paths. 

To further investigate the quality of lower and upper bounds 
and the impact of different values of Njj on the performance of 
successive convex relaxation technique, we obtain the lower 
and the upper bounds, their gaps defined as the difference 
between upper and lower bounds, and the number of iterations 
for upper bounds for the value of N from 2 to 50. Figs. 
and 1^ show the results for energy cost minimization and PAR 
minimization, respectively. 

From the figures, we observe that the overall trend in the 
results from the minimization of energy cost and PAR are 
similar to each other, except the relatively larger gaps in the 
PAR minimization that we already discussed with the results 
shown in Fig. For both minimization problems, the upper 
bounds are quite close to the lower bound for a broader range 
of values of N, and the curves for lower bounds with different 
values of Nu are hardly distinguishable. The gaps shown in 
Fig. ID (b) and Fig. (b) further confirm that the quality of 
suboptimal solutions represented by the lower bounds hardly 
depends on the values of Nd but slightly improves as N 
increases; even though both upper and lower bound values 
roughly increase as N increases, the gaps fluctuate but do not 
show any trend of increasing. 

While Fig. |^(a) and (b) and Fig.|^(a) and (b) show that Njy 
does not have any visible impact on the quality of suboptimal 
solutions. Fig. (c) and Fig. (c) show that increasing the 
value of Nu can significantly reduce the number of iterations. 
As discussed in Section |III-B| these results show that a larger 
value of Nd, in combination with a lower dropping threshold 
value (i.e., 9d = 0.1), can improve the speed of the successive 
convex relaxation while maintaining the quality of resulting 


suboptimal solutions. 

V. Conclusions 

The atomicity of appliance operation has never been given 
serious attention in scheduling of appliance energy consump¬ 
tion in autonomous DSM for residential smart grid. The 
current dominant approach based on the vector of appliance’s 
hourly energy consumption may result in several gaps in 
scheduled appliance energy consumption and have problems 
in providing enough power for appliances to carry out required 
operations, which, in most cases, are atomic. 

In this paper, therefore, we have provided a new formulation 
of appliance energy consumption scheduling based on the 
optimal routing framework, which guarantees the atomicity of 
resulting scheduled energy consumption as such without ad¬ 
ditional objective functions. Compared to the straightforward 
problem formulation based on the vector of possible starting 
times of an appliance operation, the optimal-routing-based 
formulation provides a Boolean-convex problem for a convex 
objective function, which is amenable to the successive convex 
relaxation technique where we can apply well-known interior- 
point methods for the efficient solution of relaxed convex 
optimization problem. Unlike approaches based on heuristics 
like genetic algorithms | [T9) , convex relaxation enables us to 
carry out systematic analysis of the original problem with both 
upper and lower bounds. 

The numerical results for the cost and PAR minimization 
problems demonstrate that the proposed successive convex 
relaxation technique can provide tight upper and lower bounds 
and, therefore, suboptimal solutions very close to optimal 
solutions, often identical to true optimal values, in an ef¬ 
ficient way. The results also show that, using two control 
parameters, i.e., Nd for the maximum number of fractional¬ 
valued elements that can be dropped per iteration and Od for a 
dropping threshold, we can strike the right balance between the 
quality of suboptimal solutions and the number of iterations 
to obtain them in applying the successive convex relaxation 
technique. Considering that the original problem of atomic 
scheduling is a very difficult combinatorial problem to solve 
using direct enumeration due to its huge size of the feasible set, 
the proposed successive convex relaxation technique makes 
it practical to implement the atomic scheduling of appliance 
energy consumption for autonomous DSM in residential smart 
grid. 

Note that our major focus in this paper is on the formu¬ 
lation of atomic scheduling problem and efficient solution 
techniques based on convex relaxation with adaptive dropping 
of fractional-valued elements. The extension to distributed 
atomic energy consumption scheduling and advanced tech¬ 
niques refining the feasible region at each iterative step to 
reduce the total number of iterations are interesting topics for 
further study. 
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Fig. 5. Atomic energy consumption scheduling for cost minimization: (a) Upper (UB) and lower bounds (LB); (b) gaps; (c) number of iterations. 





—• UB{1)-LB 
UB{2)-LB 
--UB{5)-LB 
- - UB{10)-LB 



UB{2) 




/■! 

-- UE 

?{in) 



/'' 


















2 10 20 30 40 50 

N 


(a) 


(b) 


(c) 


Fig. 6. Atomic energy consumption scheduling for PAR minimization: (a) Upper (UB) and lower bounds (LB); (b) gaps; (c) number of iterations. 
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